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ABSTRACT 

This  paper  presents ?an  analytic  method  for 
studying  the  mechanisms  that  lead  to  the 
detonation  of  cased  munitions  filled  with 
molten  explosives  when  Impacted  by  large 
concrete  fragments,  representing  failed  wall 
sections.  This  Is  a single  degree  of  freedom 
dynamic  structural  analyst s.  method. which 
makes  use  of  predetermined,  nonlinear  load- 
deflection  characteristics  of  the  shell 
casing;  In  the. analysis  described,  the 
method  was  applied  tb  predict  the  pressure- 
time history  ,1n  the  molten  explosive  when 
subjected  to  Impact.  A reasonable  comparison 
between  analytic  and  experimental  results  was 
obtained. 


INTRODUCTION 

At  one  state  of  the  munition  production  pro- 
cess, molten  explosives  are  poured  Into  empty  shell 
casings.  At  this,  and  most  stages  of  the  producr 
tlon  process  there  Is  a finite  probability  that  an 
accidental  explosion  may  occur.  To  limit  the 
propagation  of  an  accidental  explosion,  the  pouring 
area  Is  subdivided  Into  smaller  areas  separated  by 
reinforced  concrete  walls.  This  does  not  com- 
pletely eliminate  the  danger  of  propagation 
because  art  explosion  In  a cubicle  can  cause  some 
breakup  of  dividing  walls,  producing  energized 
fragments  which  may  Impact  shell  casings  In 
neighboring  cubicles  thus  possibly  producing 
additional  detonations.  Whether  or  not  an  explo- 
sion is  produced  by  such  Impact  depends  on  the 
dynamic  pressures  produced  within  the  molten 
explosive,  which  In  turn  depends  on  the  mass  of  the 
impacting  fragment,  its  impact  velocity,  point  and 
angle  of  the  impact,  etc.  The  size  of  fragments 
produced  by  a separation  wall  depends  on  the 
physical  characteristics  of  the  wall,  1;e.  rein- 
forcement details,  concrete  strength,  aggregate 
size,  thickness,  span,  support  conditions,  etc. 

It  also  depends  on  the  intensity  and  distribution 
of  the  blast  load  produced  by  the  donor  charge  on 
the  dividing  wall.  . 

This  paper  describes  a project  (Ref.  1)  whose 
objective. was  to  study  the  mechanisms  controlling 


the  Impact  sensitivity  of  shell  casings  filled 
with  molten  explosives  when  impacted  by  secondary 
fragments.  (It  was  determined  In  earlier  experi- 
mental investigations  that  a shell  in  the  molten 
or  just-poured  state  was: more  sensitive  to 
Initiation  than  when  the  explosive  was  cool  and 
hardened).  The  project  was  divided  Into  two 
parts. 

The  first  part  was  an  analytic  effort  whose 
purpose  was  to  predict  the  pressure  buildup  within 
the  molten  explosive  when  the  casing  and  the  con- 
tained explosive  Interact  with  a secondary 
fraoment.  It  Is  hypothesized  that  pressure  build- 
up (or  the  rate  of  pressure  buildup)  and  the 
resulting  rise  In  temperature  is  a good  measure 
of  explosive  detonation  sensitivity. 

The  second  part  was  an  experimental  effort 
whose  aim  was  to  ascertain  the  credibility  and 
accuracy  of  the  analysis.  Instrumented,  full 
scale  experiments  were  performed  to  obtain 
quantitative  data  for  comparison  with  analytic 
results.  The  experiments  also  provided  an  In- 
sight Into  what  was  occurring  during  the  impact 
process  and  thus  helped  guide  the  analysis. 

ANALYSIS  OF  SHELL  EXPLOSIONS  DUE  TO 
IMPACT  BY  LARGE  CONCRETE  FRAGMENTS 

Figure  I Is  a schematic  showing  the  shell - 
fragment  configuration  analyzed.  The  shell 
casing,  open  at  its  apex  and  completely  filled 
with  molten  explosive,  Is  stationary  on  the 
ground  plane.  At  the  time  of  Impact  the  explosive 
Is  at  a sufficiently  high  temperature  to  be  in  a 
liquid  state.  The  fragment,  a concrete  cylinder, 
is  moving  in  a horizontal  direction  with  velocity, 
V.  The.  conditions  are  assumed  to  be  such  that 
during  Impact  the  shell  casing  will  be  plastically 
deformed  and  will  experience  an  acceleration.  The 
plastic  deformation  of  the  casing  will  cause  the 
molten  explosive  to  be  pressurized  and  to  be 
forced  to  flow  up  and  through  the  filling  orifice 
at  the  apex  of  the  casing.  This  pressure  buildup 
during  the  time  of  impact  and  the  associated 
temperature  rise,  is  hypothesized  to  be  a measure 
for  determining  if  the  explosive  will  detonate.. 

The  analytic  procedure  used  to  obtain  the  pressure- 
time relationship  in  the  explosive  during  impact 
is  described  next. 


44- 


i 


I 

) 


} 


{ 

! 


\ 


A 


FIG.  I.-  She! i -Fragment  Configuration 


Analytic  SolutionRegulreroents 

The  analytic  solution  Is  a single-degree  of 
freedom  dynamic. analysis  which  requires  certain 
postulated  explosive,  shell  casing,  and  frapent 
Interaction  characteristics  which  are  the 
following:  . <■ 

1.  Force-deflection  characteristics  (resistance 
function)  of  the  shell  casing 

2.  Volume  versus  deflection  characteristic  of 
the  shell  casing 

3.  Orifice  area  versus  deflection  characteristics 

■ For  the  purpose  of  determining  the  forces 
deflection  characteristics,  the  shell  casing  was 
assumed  to  be  restrained  In  the  horizontal  direc- 
tion*ohly, all  along  the  line  formed  by  the  Inter- 
section of  the  midplane  of  the  shell  surface  and 
the- vertical  plane  whose  edge  view  is  Indicated  by 
lineA-A  In  Fig.  I.'1  Making  use  of  the  .shell 
symmetry,  one  quarter  of  the  shell  surface  Is 
modeled  for  finite  element  analysis  as  shown  in 
Fig.  2.  Thickness  is  varied  along  the  height  In 
discrete  Intervals.  AHSYS  (Ref,  2)  finite  element 
computer  program  which  meets  the  requirements  of 
large  deflection  and  plastic  deformation  of  the 
casing  material,  was  used  in  the  analysis.  Force- 
deflection  characteristics  were  obtained; by 
Imposing  deflections  and  then  computing  the  corre- 
sponding forces  required  to  produce  them.  Deflec- 
tions were  imposed  on  the  shell  casing  at  nodes 
representing  ah  area  approximately  equal  to  the 
contact  area  between  the  casing  and  the  impacting 
concrete  cylinder.  Load  deflection  characteristics 
were  determined  for  three  pressures  applied  to  the 
interior  of  the  shell.  Results  are  shown  in. Fig. 

3.  Along  with  force-deflection  characteristics  of 
the  shell  casing,  volume-deflection  and  orifice 
area-deflection  characteristics-  were  obtained 
from  this  analysis;  These  results  are  shown  in 
Fig.  4 and  Fig,  5- respectively-.  It  will  be  noted 
that  both  of  these  characteristics  are  Independent 


of  Internal  pressure.  Also,  at  least  for  this 
shell  casing;  there' was  essentially  no  change  in 
the  orifice  area  with  change  in  deflection,  see 
Fig.  5. 
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FIG.  2. -Finite  Element  Mesh  of  One  Quarter  of  the 
Shell 
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FIG.  3. -Force  Versus  Deflection  Characteristics 
of  the  Shell  Casing  for  Indicated 
Internal  Pressures 
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FIG.  4. -Volume  Versus  Deflection  Characteristic 
of  the,  Shell  Casing 


In  the  Initial  analysis  It,  was  .assumed  that  the 
molten  explosive  was  Incompressible.  Analytic 
results  using  this  assumption  showed  large  rapid 
changes  In  the  liquid  explosive  pressure;  suggest- 
ing compressibility.  The  analysis  was  modified  to 
take  Into  .account  fluid  compressibility.  In  the 
modlflediahalysls  the  characteristics  shown  In 
Fig.  6 are  still  used,  however  the;  corresponding 
Internal  pressure  Unmodified  .on 'the  basis  of  the 
bulk  modulus  of  the  explosive,  i;e. 

B,,p0  ; (2) 

where  B * the  bulk  modulus  taken  as  580*000  psl 
p ■ the  weight  density,  of  the«molten 
explosive 

ap  * change  In  Internal  pressure 
Analytic.  Solution  Procedure 


FIG-.  5. -Orifice  Area  Versus  Deflection 

Characteristic  of  the  Shell  Casing 

Next,  we  were  concerned  with  the  flow 
characteristics  of  the  mpltep?explos1ve  through 
the  filling  orifice  of  the  casing.  The  flow  Is 
dependent- on 'the -orifice  area,  the  fluid  ^pressure, 
viscosity  and  the  orifice  coefficient.  The 
following  formula  (Ref.  3)  was  used  to  relate  the 
flow  rate,  Q to  the  Internal  pressure,  p of  the 
liquid  explosive  and  other  parameters. 


where  Cc  = the  orifice  coefficient,  taken  as  0.64 


A * the  orifice  area 
g = acceleration  due  to  gravity 
Y * weight  density  of  the,  molten  explosive, 
taken  ai  a'pproximately  85  lb/cu  ft 

Flow  character! sties, are  shown  in  Fig.  6 for  three 
orifice  areas. 


As  lndlcat'ed. earlier,  the  analytic  solution  Is 
a single-degree  of  fr,eedom<analysis.  Referring 
to  Fig,  1,  assuming  that  the  concrete  cylinder  and 
the  shell  move  -In  a straight  line*,  the  equations 
of  motion  and  Initial  conditions  are -given  as 
follows; 

M.  V-  V-F  (3) 

Ms  k;  - F (4) 

where  Mc  and  Ms  represent  the  mass  of  the  cylinder 
and  the  shell  casing  respectively..  F is  the  inter- 
act1on--force> between  the  shell  and  the.  concrete 
cylinder,  and  x&  ahd  xV  are  the  center  of  gravity 
accelerations  of  the  concrete  cylinder  and  the 
shell  respectively. 

The  Initial  conditions  for  the  concrete 
cylinder  and  the  shell-  are  (at  t*0.  Impact  Is 
Initiated)  x<-;  the  Initial:  veloclty-of  the, cylin- 
der Is  equal  to  V;  xs»  the  Initial:  velocity  of  the 
shell  is  equal  to  zero.  The  initial- displacements 
xc  and  xs,-of  the  cylinderand  the  shell  are  .'equal 
td -zero.  Euler  equations  were  used  in  the  Inte- 
gration? process*  1»e. 
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X (t  + At)  = *x(t)  + *x(,t)At  (5)= 

x (t  + At)  = x(t)  + x(t)At  (6): 


The  solutl on  procedure  is  detailed  below. 

SolutionProcedure 

1.  Read-data  describing  the  system 

2..  Set  time  to  zero  (t=0)  and  displacements  and 
ve16c1ti.es  to  their  Initial  values.  Set 
initial  internal  pressure  to  zev-b,  and  the 
orifice  area  and  the  casing  volume  to  the 
Initial  values.* 

3.  Printout  time  and  pressure 

4.  If  time  exceeds  maximum  value,  stop., 

5.  From  force'-def  lection-pressure  curves  (Fig.  3) 
determine  the  Interaction  force,  F.  Note 
that  the  casing  deflection  <5,  Is  the  differ- 
ence in  the  cylinder  and  shell  motions 
(displacements),  i.e.  S * xc-xs. 

6.  Compute  cylinder  and  shell  accelerations,  see 
equations  (3)  and  (4). 

T.  Use  the  Euler  Integration  formulas, 
expressions;  (5)  and  (6),  to  determine 
velocities  and  displacements  of  the  cylinder 
and  the  casing  at  time  t + At.  Update  time 
to  t + At. 

8.  Determine  shell  casing  volume  using  volume- 
deflection-pressure  curve  (Fig.  4).  Compute 
volume  rate  of  flow  from  difference  In- volume 
from,  previous  time  step  .divided  by  time  step, 
At. 

9.  Determine  orifice  area  from  orifice  area- 
clef  lection-pressure  curves  (Fin,  5). 

10.  Determine- new  Internal  pressure  from 
pressure-volume- flow  rate-area  curves 
(Fig.  6). 


Steps  11  through  13  represent  an 
Iterative  procedure  used  to  take 
Into  account  the  compressibility 
of  the  liquid  explosive. 

ili  Determine  the  change  in  pressure  over  the 
time  step  At  and  compute  the  change  in  the 
weight  density  Of  the  liquid  explosive,  see 
Eq.  (2), 

12.  Determine  the  new  weight  density, 

13.  Determine  the  change  In  volume  duetto  the 
change  In  the  weight  density.  Compute?  the 
volume  rate  of  flow. 

Continue  the  iteration  until  tne 
change  in  volume  is  smaller  than 
a preassigned  value, 

14.  Return  to  step  3 and  continue. 

EXPERIMENTAL  EFFORT 

The  purpose  of  the  experimental  effort  was  to 
measure  the  pressure-time  history  within  the 
molten  explosive  when  the  shell  containing  it  is 
impacted  by  a concrete  fragment.  The  experi- 
mental setup  is  depicted  in  Fig,  7.  A shell 
casing,  containing  a liquid,  Glycerol,  of  the 


same  density  :as  molten  cbmposition  'B  explosive 
is  mbunted'^jn^i^y^bh' :S¥peide'stal£ -is'eS 'itatrget"' ' 
in  Fig.  7.  Ah  air  gun  fs  then  used  to  launch 
the  .concrete  fragment  at  the  shell  casing  to  impact 
at  a specified  aiming  point;  The -concrete*  fragment 
IS  In  the  slape  of  a -cylinder’.  ^Two’sizes  and 
weights  .were *used  in- this  study  (Ref.  1);  2-ft  and 
4-tf.t  long. and 'weighing  200- lbs  and  400-lhs  respec- 
tively. The  InstriiwnentatTon  c'onsi sted^of  a prelsf> 
sure  sensor  located?W1tH1h  the  liquid  for  measuring 
the  pressure-time' history  during  impact.  Photo- 
graphic coverage  was  provided  to  measure  the  impact 
velocity  of  the  concrete  cylinder... 


FIGt  7.-Secondary;Fragments  Impact  Test  Site 


Selected  Results 

Figure  8 shows  an  experimentally  determined  and 
the  corresponding  analytic  pressure-time  history. 
This  particular  experiment  dealt  witlua  4,2-1n*. 
mortar  sly?1!  filled  with  Glycerol  and  water  at 
ambient  |e  >perature.  A concrete* cylinder,  simulat- 
ing a Wat  1 fragment  and. weighing  200  lbs. , impacted 
the  shelr  -t  272  ft/sec  at  a point  4;25  in  below 
the  filVimrorifice.  In  terms  of  the  shape  of  the 
pressure-time  history*  duration  and  peak  pressure* 
the  comparison  between  experimental  and  analytic 
results  appears  to  be  favorable.  Similar  compari- 
sons vieM?  obtained  for  other  experiments  conducted 
in  the  mirse  of  the.  study. 

CONCLUSIONS  AND  RECOMMENDATIONS 

A simple  analytic  method  was  formulated  for 
predicting:  pressure-time- hi  stories  in  shell- 
casings  filled  with  molten  explosives  when.  Im- 
pacted by  secondary  fragments.  Analytic  results 
compared  favorably  with  experiments. 

Additional  analyses  are  required  to  study  the 
sensitivity  of  results  for  a larger* set  of  experi- 
ments with  the  objective  of  improving  the  accuracy 
of  the  predictive  method.  To  date  four  other  ? 
experimental  programs  have  been  conducted  on  the 
sensitivity  to  Impact,  by  large  concrete  fragments, 
of  a variety  of  molten  and  ambient  temperature 
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explosive-filled, shells.  Fotvthese. previous 
experimental,  programs  apaly  ti cal  >predi ctl ons  of 
tlnie-  cM.tv^sji(^1d^|)erfom^i.  v[it^.  a larger 
number  of  curves,.. andhehce  larger  number  of  test, 
icpnditlonsi.  one  would'be  able  to  distinguish 
betwepvpressure-tline  cpndl tioris  for  explosion  :and 
no,;  reaction,  The  scope  of  thls  method  needs  to  be 
further  expanded  to  study  whether  the  peak>, pressure 
in  'the  liquid  (explosive,  rate  of  pressure  rise, 
etc.  are  parameters  v*1  ch; individually,  or  in  com- 
bi nation  will  predi ct  the  onset  of  detonation. 


FIG,  8. -Pressure  Versus  Time  Graph  for 

Experiment  4 (Concrete  Projectile  200  lb, 
at  272  ft/sec  Velocity) 
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